Muon Spin Relaxation Investigation of the Spin Dynamics of Geometrically 
Frustrated Antiferromagnets Y 2 Mo 2 07 and Tb 2 Mo 2 07 
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The spin dynamics of geometrically frustrated pyrochlore an- 
tiferromagnets Y2M02O7 and TD2M02O7 have been investi- 
gated using muon spin relaxation. A dramatic slowing down 
of the moment fluctuations occurs as one approaches the spin 
freezing temperatures (TV =22 K and 25 K respectively) from 
above. Below Tj? there is a disordered magnetic state similar 
to that found in a spin glass but with a residual muon spin 
relaxation rate at low temperatures. These results show that 
there is a large density of states for magnetic excitations in 
these systems near zero energy. 
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Antiferromagnets which are frustrated or diluted can ex- 
hibit novel electronic and magnetic behaviour. Recently, 
there has been considerable interest in the behaviour of 
systems where the natural antiferromagnetic coupling be- 
tween ions is frustrated by the geometry of the lattice. In 
two dimensions, Heisenberg spins on triangular and cor- 
ner sharing triangular (kagome) lattices are simple ex- 
amples of geometric frustration, while in three dimen- 
sions, the most well studied systems have a pyrochlore 
structure, in which the magnetic ions occupy a lattice of 
corner sharing tetrahedra. A system of Heisenberg spins 
interacting via nearest-neighbor antiferromagnetic cou- 
plings on the pyrochlore lattice displays a classical ground 
state with macroscopic degeneracy, since the lowest en- 
ergy spin configuration requires only that £ 4 =1 Si = 
for each tetrahedron. This feature led Villain to argue 
that these systems remain in a cooperative paramagnetic 
state with only short range spin-spin correlations for all 
T>0 |lj and this has been confirmed by Monte Carlo 
simulations Possibly, the most interesting feature 
of the ground state of pyrochlore || and kagome 
lattice antiferromagnets is the prediction of a dispersion- 
less spin- wave branch ("zero modes"). These zero modes 
manifestly affect the thermodynamics of these classical 
systems as demonstrated by Monte Carlo simulations, 
where the low temperature specific heat, C v , falls below 



the classical value ks expected from equipartition of en- 
ergy HHH. Ah* ' again because of these zero modes, the 
kagome and pyrochlore antiferromagnets display large 
spin fluctuations down to T= + @. However, further 
nearest neighbour exchange || , magnetic anisotropy and 
fluctuations ||] may lift this classical ground state degen- 
eracy. 

A wide variety of interesting magnetic behaviour has 
been observed in real systems. Neutron scattering re- 
sults on FeF3 || show a transition to a non-collincar 
long range ordered (LRO) state, in which the spins on 
a tetrahedron point away from the center. However, 
a large number of oxide pyrochlores do not show Neel 
LRO. Bulk magnetic susceptibility measurements |lC| ] 
on the pyrochlore Y2M02O7 show strong irreversible be- 
haviour below Tp=22 K, characteristic of spin glass or- 
dering, even though the level of disorder is immeasurably 
small. Recent measurements of the dc magnetization of 
Y2M02O7 show a divergent non-linear susceptibility at 
Tk 22 K pi , which is a signature of a true thermo- 
dynamic spin glass phase transition. Inelastic neutron 
scattering data |Q on Tb2Mo207 confirm there is rapid 
slowing down of the Tb spins as one approaches T^=25 
K from above and the absence of LRO below Tp. The 
observed strong diffuse scattering in Tb2Mo207 indicates 
the presence of short range correlations between the mo- 
ments, which are frozen on a time scale of about 10~ n 
s. 

In this letter we report an investigation of the low tem- 
perature magnetic properties of pyrochlores Y2M02O7 
and TD2M02O7 using the technique of muon spin rota- 
tion/relaxation (//SR ), which is sensitive to spin fluctu- 
ation rates in the range 10 4 — 10 11 s~ 4 below that 
detectable with neutron scattering. We find that, de- 
spite its nominally disorder free structure, the magnetic 
behaviour in Y2M02O7 is close to that observed in con- 
ventional random spin glasses. Specifically, a large static 
internal magnetic field with a very broad distribution de- 
velops below Tp, such that no coherent muon spin pre- 
cession is observed. At the same time, the muon spin 
relaxation rate 1/T\ decreases according to a power law 
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with decreasing temperature. A similar magnetic transi- 
tion occurs in TD2M02O7. The most remarkable feature 
in the data is the presence of a sizeable residual spin re- 
laxation rate at low temperatures, which is not evident 
from previous data on conventional metallic spin glasses 
like CuMn @, AuMn |l5j and amorphous-FeMn |6|. 
This is direct evidence for a larger density of magnetic 
excitations near zero energy than in conventional random 
spin glasses. 

Details on the preparation of the Y2M02O7 and 
Tb2Mo2C>7 samples are given elsewhere |10 . Pyrochlores 



ically equivalent muon site i can be described by a single 
exponential e~ Xit with a relaxation rate [0 : 



crystallize with an fee structure containing eight formula 
units per conventional unit cell and space group Fd3m. 
The ions on the 16d site form a network of corner sharing 
tetrahedra; the 16c sites constitute an identical sublat- 
tice, displaced by Mo 4+ ions occupy the 16c 

site, Y 3+ or Tb 3+ ions the 16d site. The Tb 3+ ion has 
a large magnetic moment of ~ 9/ib > roughly nine times 
larger than that of Mo 4+ , whereas Y 3+ is diamagnetic. 
Y2M02O7 and Tb2Mo207 are semiconductors with small 
band gaps of 0.013 and 0.007 eV respectively @- The 
samples in this study were characterized by magnetic sus- 
ceptibility and x-ray diffraction. Sharp irreversibilities in 
the magnetization were observed at spin freezing temper- 
atures of 22 K and 25 K respectively, consistent with that 
seen in other highly stoichiometric samples of Y2M02O7 
and TD2M02O7. ^From Rietveld profile refinements of 
neutron diffraction measurements one can say the 
concentration of oxygen vacancies, likely the main source 
of crystalline disorder in these materials, is below the de- 
tectable limit of 1%. Since the ionic radii of Y 3+ and 
Mo 4+ , as well as that of Tb 3+ and Mo 4+ are very differ- 
ent, there should be no admixing between the 16c and 
16d cations. This is confirmed by analysis of X-ray data. 

fiSR measurements were made at TRIUMF in a 4 He 
gas flow cryostat for temperatures above 2 K and in an 
Oxford Intruments Model 400 top loading dilution refrig- 
erator (DR) for lower temperatures. For the DR mea- 
surements the pressed polycrystalline pellets were var- 
nished onto an Ag plate and covered in thin Ag foil, which 
was bolted to the cold finger. In a /iSR experiment the 
observed quantity is the time evolution of the muon spin 
polarization, which depends on the distribution of inter- 
nal magnetic fields and their temporal fluctuations. In 
a longitudinal field (LF) geometry an external magnetic 
field is directed along the initial polarization direction. 
The present measurements were made in a small longi- 
tudinal field to quench any spin relaxation from static 
nuclear dipolar fields in the sample holder. Further de- 
tails on the ^iSR technique may be found in Ref. |l3|] . 

Fig. 1 shows several typical /iSR spectra in Y2M02O7. 
Above Tp=22 K the observed spin relaxation is at- 
tributed to rapid fluctuations of the internal magnetic 
field due to Mo 4+ moments in the paramagnetic phase. 
When the fluctuation rate v A (defined below), the 
relaxation function [see P z (t) in Fig. 1] for each magnet- 
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where Aj = 7 M i?i is the gyromagnctic ratio of the muon 
(2ir x 135.54(10 6 rads s _1 tesla -1 )) times the rms inter- 
nal magnetic field Bi at site i. i>i is the fluctuation rate 
of the internal field and vl — 7 M -B C xt is the Larmor fre- 
quency of the muon in the external magnetic field. Note 
that \i is only weakly dependent on the applied field pro- 
vided Ui ^S> Vl\ this is consistent with the absence of any 
field dependence observed in the spectra for T above Tp. 
Fig. 2 shows the average muon spin relaxation rate in 
Y2M02O7 obtained from fits to a single exponential re- 
laxation function P z (t) ~ e~*/ Tl over a restricted time 
interval of 0.05 to 6 fj,s, where A = 1/Ti. In the param- 
agnetic phase one may use Eqn. 1 to estimate the aver- 
age fluctuation rate of the moments. For example with 
Bi = 0.066 T (see below) one obtains fluctuation rates 
shown in the inset of Fig. 2. Note the sharp rise in the 
average 1/Ti and corresponding decrease in the Mo 4+ 
fluctuation rate as one approaches Ti?=22 K. 

Just above Tp, P z (i) deviates somewhat from a single 
exponential (see for example T=27.5 K spectrum in Fig. 
1) and is better described by a stretched exponential of 
the form e~( At ) , with (3 near 0.4. Similar behaviour has 
recently been observed in other dense spin glasses AgMn 
and AuFe @. 

The muon spin polarization function below Tp (see 
inset in Fig. 1) is characterized by rapid depolarization 
of 2/3 of the initial polarization, followed by slow relax- 
ation of the remaining 1/3 component. This is a char- 
acteristic signature of a highly disordered magnetic state 
in which the moments are quasi-static on the timescale 
of the muon lifetime. For example, the muon polariza- 
tion function for a single magnetic site with a Gaussian 
distribution of static internal fields is given by the Kubo- 
Toyabe function |TJ]: 
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The curve in the inset of Fig. 1 shows a fit of the early 
time data at 2.5 K to Eqn. 2, modified slightly to in- 
clude the small external field of 0.02 T. The best fit gives 
a value A/7 M = 0.066(3) T, which corresponds to an av- 
erage field strength n /87ttA/7 ai = 0.105(5) T. Note how- 
ever that the dip in P z (t) at 0.032 fis is not as deep as 
predicted by the modified Eqn. 2, indicating the distri- 
bution of internal fields is more complicated than a single 
Gaussian. One can generalize Eqn. 2 to include a fluc- 
tuating component to the internal field, which results in 
relaxation of the 1/3 tail seen in Fig. 1. 

Muon spin relaxation results from the exchange of en- 
ergy with magnetic excitations. A first order process, in 
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which the muon absorbs or creates an excitation with an 
energy equal to the muon Zeeman energy, is normally 
suppressed in conventional systems with LRO, where the 
density of states p(E) — » as E — > 0, since it requires 
excitations near zero energy. In a second order (Raman 
magnon scattering) process involving inelastic scattering 
of an excitation, application of Fermi's Golden rule gives: 
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where the muon Zeeman energy has been neglected and 
M(E) is the matrix element for inelastic scattering of an 
excitation of energy E causing a muon spin flip. In a spin 
glass, n(E /ksT) is the probability distribution (assumed 
to be Bose) for "intravalley" excitations, ie. spin excita- 
tions within one of the macroscopic number of metastable 
states or valleys. ^From Eqn. 3, the temperature depen- 
dent behaviour of 1/Ti is primarily determined by the 
energy dependence of p(E)M(E). The low temperature 
linear specific heat observed in Y2M02O7 suggests 
p(E) is flat or at least weakly dependent on energy. If 
p(E) and M(E) have power law dependences with pow- 
ers I and m respectively, then Eqn. 3 implies that 1/Ti 
varies as 7 l2 ('+™ l )+ 1 below TV- In other words, below 
Tf , 1/T\ decreases gradually as the magnetic excitations 
freeze out. The curve in Fig 2 shows the best fit of the 
data below 12 K to a simple power law form A = Xo+AT n 
with exponent n = 2.1(3). This power law behaviour and 
the small value of n indicate that M(E)p(E) in Eqn. 3 
has a very weak energy dependence (ie: l + m is less than 
1). Intervalley transitions, involving reorientations of fi- 
nite sized spin clusters, are thought to be important only 
in the mK range fl9|l , where 1 jT\ is independent of tem- 
perature in this sample. There is a small residual relax- 
ation rate ( Ao =0.02 /is -1 ) at the lowest temperatures, 
which implies there is a non-zero density of excitations 
close to zero energy. Such relaxation is just above the 
resolution limit of the /iSR technique. 

A similar spin freezing transition is observed in 
Tb 2 Mo207, but the residual 1/Ti at low temperatures 
is much larger. Fig. 3 shows the muon spin relaxation 
rate measured in a small longitudinal field of 5 mT. As 
in Y2M02O7, a critical slowing down of the moment fluc- 
tuations occurs as one approaches Tf =25K from above. 
Using a value for Bi = 0.7T (see below) we obtain the 
spin fluctuation rates (u) above Tf shown in the inset 
of Fig. 3. For comparison, we include some of the cor- 
responding Tb 3+ spin fluctuation rates determined from 
inelastic neutron scattering p2[ . Considering the fluctu- 
ation rates measured by neutron scattering are at the 
lower experimental limit and the systematic errors in 
both measurements, the agreement is reasonable. ^From 
this we can conclude that both techniques are sensitive 
to the same quantity in this sample, ie: the Tb 3+ mo- 



ment fluctuation rates. The fact that Tf is about the 
same in Tb2Mo20y and Y 2 Mo 2 07 supports the proposal 
that the spin freezing temperature in Tb2Mo207 is de- 
termined mainly by the Mo 4+ ions, which provide an 
effective coupling between the larger but more localized 
rare earth Tb 3+ moments. 

Fig 3 shows that 1/Ti initially begins to decrease as T 
falls below T p but recovers below 1 K and stays constant 
at a relatively large value of 5 ps -1 . The initial ampli- 
tude of the relaxing 1 /3 component increases as the ratio 
between the external magnetic field and internal static 
field. This dependence was used to estimate the mag- 
nitude of the static component of the internal magnetic 
field Bi = 0.70(6) T which is about an order of magni- 
tude larger than in Y2M02O7 |2(J, as expected from the 
ratio of Tb 3+ and Mo 4+ magnetic moments. This con- 
firms that the Tb 3+ moments are involved in the 25 K 
freezing transition. Note the ratio of residual relaxation 
rates in Y2M02O7 and Tb2Mo207 is roughly equal to 
the ratio of the square of the respective internal fields. 
The large residual 1/Ti in Tb2Mo2C>7 establishes there 
is a non-zero density of low energy excitations, which 
cause relaxation either by a first or second order process. 
Computer simulations by Ching et al. pl[ | on insulating 
Heisenberg spin glasses Eu^Sri^^S (x = 0.54 and 0.40) 
have indicated the density of states p(E) may be peaked 
at low energies and p(0) finite. 

We emphasize here that we find convincing evidence 
for a limiting temperature independent 1/Ti in Y2M02O7 
and Tb2Mo207 only in the temperature range T/Tf < 
0.05. Previous pSK experiments |l4|-|l6[] found a strong 
temperature dependence of 1/Ti in the temperature 
range T/Tp e [0.1 — 1.0], with no sign that 1/Ti was ap- 
proaching a limiting and temperature independent value 
linix^o [1 /T\ (T)] above the experimental ^SR resolution 
limit and in any case, did not probe the temperature 
range T/Tf < 0.1. It is interesting to note that other 
spin glasses like Cdi-^Mn^Te (0.27< x < 0.65) |§§] and 
Lai.g4Sro.o6Cu04 p3|] show indications of low tempera- 
ture spin dynamics but again, these insulating Heisen- 
berg spin glasses have not been studied in the impor- 
tant region below O.lTp . The geometrically frustrated 
kagome lattice system SrCrsGa40i9 has also recently 
been studied using /iSR. Dynamics spin fluctuations are 
observed without static freezing, even at lOOmK, well be- 
low Tf=3.5 K [p4j| . There is however some controversy 
over SrCrgGa40ig pE,M, as it has been suggested that 
this material does not show a thermodynamic freezing 
transition at Tf [^6). In this case, one would expect to 
find spin dynamics persisting down to zero temperature. 
This is not the case for Y2M02O7 where we have strong 
evidence for a collective freezing transition at Tf as seen 
in the critical slowing down seen in /iSR and the diver- 
gent nonlinear susceptibility pT| . 

In conclusion, despite the nominal absence of disor- 
der, the freezing process in Y2Mo 2 07 and Tb 2 Mo207 
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appears similar to that expected for a dense spin glass. 
In particular we observe a critical slowing down of the 
spin fluctuations and non exponential muon spin relax- 
ation near TV, while below Tp there is evidence for a 
highly disordered magnetic structure. The most strik- 
ing feature in both systems is the presence of a residual, 
temperature independent spin relaxation which persists 
down to very low temperatures. This shows there is an 
appreciable density of states for low energy magnetic ex- 
citations which is much larger in these systems than in 
conventional randomly frustrated spin glasses. It is pos- 
sible that the residual low temperature dynamics in these 
systems are "remnants" of the zero-modes predicted the- 
oretically for nearest-neighbor Heisenberg spins on a py- 
rochlore lattice 

This research has been partially funded by the NSERC 
of Canada under the NSERC Collaborative Research 
Grant Geometrically- Frustrated Magnetic Materials. 



[25] A. P. Ramirez, G. P. Espinosa, A. S. Cooper, Phys. Rev 

Lett. 17, 2070 (1990). 
[26] B. Martinez, A. Labarta, R. Rodriquez-Sola, X. 

Obradors, Phys. Rev B 50, 15779 (1994). 

FIG. 1. The muon spin relaxation function, P z (t) at var- 
ious temperatures in Y2M02O7. The inset shows the early 
time behaviour at T — 2.5 K. 

FIG. 2. The muon spin relaxation rate 1/Ti vs tempera- 
ture for Y2M02O7 in a small applied field of 0.02 T. The solid 
line is the best fit to the data assuming a power law functional 
form. The Mo 4+ spin fluctuation rate vs temperature above 
Tf is shown in the inset. 

FIG. 3. Dynamical muon spin relaxation rate 1/Ti vs 
temperature for Tb2Mo207 in an applied field of 5 mT. The 
inset shows the Tb 3+ moment fluctuation rate vs temperature 
above Tf- The inverted triangles indicate neutron scattering 
data [8]. 



[1] J. Villain, Z. Physik B 33, 31 (1979). 

[2] J. N. Reimers, Phys. Rev. B 45, 7287 (1992). 

[3] J. N. Reimers, A. J. Berlinsky, A. -C. Shi, Phys. Rev. B 

43, 865 (1991). 
[4] A.B. Harris, C. Kallin and A.J. Berlinsky, Phys. Rev. B 

45, 7536 (1992). 
[5] J. T. Chalker, P. C. W. Holdsworth and E. F. Shender, 

Phys. Rev. Lett. 68, 855 (1992). 
[6] J. N. Reimers and A. J. Berlinsky, Phys. Rev. B 48, 9539 

(1991). 

[7] A. Keren, Phys. Rev. Lett. 72, 3254 (1994). 

[8] S. T. Bramwell, M. J. P. Gingras, J. N. Reimers, J. Appl. 

Phys. 75, 5523 (1994). 
[9] J. N. Reimers, J. E. Greedan, M. Bjorgvinsson, Phys. 

Rev. B 45, 7295 (1992). 
[10] N. P. Raju, E. Gmelin, R. K. Kremer, Phys. Rev. B 46, 

5405 (1992). 

[11] M. J. P. Gingras et al, J. Appl. Phys. 79, 6170 (1996). 

[12] B. D. Gaulin et al, Phys. Rev. Lett. 69, 3244 (1992). 

[13] S. F. J. Cox, J. Phys. C 20, 3187 (1987). 

[14] Y. J. Uemura et al, Phys. Rev. B 31, 546 (1985). 

[15] H. Pinkvos et al, Phys Rev. B 41, 590 (1990). 

[16] M. J. P. Gingras et al, TRIUMF preprint (1995). 

[17] M. A. Subramanian, G. Aravamudan, G. V. Subba Rao, 

Mat. Res. Bull. 15, 1401 (1980). 
[18] J. A. Campbell et al, Phys. Rev. Lett. 72, 1291 (1994). 
[19] W. Y. Ching, D. L. Huber, Phys. Rev. B 34, 1960 (1986). 
[20] Since the internal magnetic fields in Tb2Mo207 are much 

larger that in Y2M02O7 it was not possible to resolve the 

fast 2/3 component and thereby measure A directly. 
[21] W. Y. Ching, D. L. Huber, K. M. Leung, Phys. Rev. B 

21, 3708 (1980). 
[22] E. J. Ansaldo et al., Phys. Rev. B 38, 1183 (1988). 
[23] B. J. Sternlieb et al, Phys. Rev. B 41, 8866 (1990). 
[24] Y. J. Uemura et al, Phys. Rev. Lett. 73, 3306 (1994). 



4 



CD 



CD 
Q 
C 
CD 

Sf 



o 
o 



o 
b 

CD 



o 
o 

03 

ro 



ro 
o 



03 

o 

o 
o 



1/T 1 (/xs 



o o o 
o o o 

- M ^ 



o o o 
^ Kj ^ 



-B 



Fluctuation Rate (THz 
o 



o 
o 
o 



o 

b 
o 



o 
b 



o 




Ejjpl 



e 



o 



~i 1 1 1 1 1 1 r 



I I I I I 



J I I I I I I I I 



J I I I I I I I I 



J I I I I I I I I 



o 
o 

o 
o 



O 

o 



1/T 1 (/uLS 



o o 

o b 
ro 



o o o 

^ ho 4^ 



o o 





o 


- en 


— I 


ro 


o 


CD 




\l 


3 


o 




~o 






CD 






~^ 






Q 


ro 






ro 


d 




ro 






- cn 


CD 




















ro 






o 






on 






o 





Fluctuation Rate (THz 

° o - 

° o o 

ro 
on 



III II 

— <g — b- 


i i 






— s- □ 




-s- □ 






□ 




□ 




□ 




□ 




□ _ 



-B- 



-B- 



-B- 



B 



-B- 



■fi- 



B 



B 



B 



-S- 



-B- 



B 



B 



B 



B 



B 



B 



on 

O 



i i i i i i 1 1 i i i i i i i 1 1 i i i i i i i 1 1 i i i 



